Effect of grain boundary on the buckling of graphene nanoribbons 
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The buckling of graphene nano-ribbons containing a grain boundary are studied using atom- 
istic simulations where free and supported boundary conditions are invoked. We found that when 
graphene contains a small angle grain boundary the buckling strains are larger when the ribbons 
with free (supported) boundary condition are subjected to compressive tension parallel (perpendic- 
ular) to the grain boundary. The shape of the deformations of the buckled graphene nanoribbons 
depends on the boundary conditions and the presence of the grain boundary and the direction of 
applied in-plane compressive tension. Large angle grain boundary results in smaller buckling strains 
as compared to perfect graphene or to a small angle grain boundary. 



Large area graphene sheets have been grown on metal- 
lic foils which were found to contain grain boundaries [1[ . 
Scanning tunneling microscopy experiments were used to 
investigate tilt grain boundaries in graphite [2j. Trans- 
mission electron microscopy was able to observe indi- 
vidual dislocations in graphene J3J. The effect of grain 
boundaries on both electronic and mechanical proper- 
ties of graphene was investigated theoretically by Yazyev 
et al [4]. Depending on the grain boundary structure, 
high transparency and perfect reflection of charge car- 
riers over remarkably large energy ranges was reported 
by using first principle calculations [5]. The effects of 
grain boundaries on the buckling of nano-scale graphene 
nanoribbons (GNR) has not been investigated up to now, 
while it is important for their mechanical stability. 

Recently, we studied the effect of applied external ax- 
ial stress on the thermomechanical properties of perfect 
graphene (PG) p, 01 and GNRS containing randomly dis- 
tributed vacancies [8] . In this paper we address the effect 
of the presence of grain boundaries and in-plane bound- 
ary compressive stress applied in two different directions, 
on the buckling and the stability of GNR for the case of 
supported boundary (SBC) and free boundary conditions 
(FBC). We found that the presence of grain boundaries: 
1) alters the sine wave shape of the longitudinal defor- 
mation modes of GNR when subjected to SBC, 2) re- 
sults in a concave shape of the GNRs when subject to 
FBC, 3) when subject to compressive tension perpendic- 
ular to the grain boundary, the GNR buckles at smaller 
(larger) strains as compared to perfect graphene in case 
of FBC (SBC), 4) the buckling transition found for com- 
pressive tension parallel to the grain boundary is three 
times larger than the one for perfect graphene subjected 
to compressive tension along the zig-zag direction, and 5) 
free energy calculations reveal that a larger angle grain 
boundary in suspended graphene when it is subjected to 
compressive tension makes it unstable. 

Initially the coordinates of all atoms are put in a flat 
surface of a honey-comb lattice with nearest neighbor dis- 
tance equal to a^ = 0.142 nm. Our perfect GNR (PG) is 
a rectangular GNR with dimensions a x b = 20 x 10 nm 2 
in x and y directions with armchair (ac) and zigzag (zz) 
edges, respectively. A grain boundary is introduced as 



FIG. 1: (Color online) Two snap shots taken from the central 
portion of GNRs with grain boundaries after relaxation at 
room temperature. Two different configurations of the grain 
boundaries (LAGBI in (a) with 6 =21.8° and LAGBII in (b) 
with =32.2°, see Ref. [j) are shown by the colored pentagon 
and heptagons . 



an array of 5-7 defects, which are put in the center of the 
PG along the y-direction with angle (see Figs. [1]). As 
an example, we study two kinds of grain bound- 
aries which were named LAGBI and LAGBII by 
Yazyev et al [4]. These grain boundaries are 
typical interfaces between domains of graphene 
with different crystallographic orientation. Mu- 
tual orientations of the two crystalline domains is 
described by the misorientation angle which for 
LAGBI is =21.8° and for LAGBII is =32.2°. 
In Figs. [1] we depict two snap shots of the cen- 
tral portion of LAGBI (a) and LAGBII (b) after 
relaxation at room temperature. 

Classical atomistic molecular dynamics simulations 
(MD) are employed to simulate compressed PG, LAGBI 
and LAGBII using Brenner's bond-order potential [9| and 
temperature is controlled by a Nose'-Hoover thermostat 
at room temperature. Before starting the compression, 
the systems are equilibrated during 75 ps (150 000 time 
steps). Extra atoms were added to the bound- 
aries of the rectangular samples which are char- 
acterized by: x = ±a/2 ± 2l, y = ±b/2 ± 2l. 
The compression and boundary conditions (FBC 
and SBC) are applied on these extra atoms which 
are outside the main systems. Compressing di- 
rection is defined by the angle 'a\ For example 
a = (tt/2) implies that compression is applied 
in the x-direction (^/-direction) so that the right 



(up) longitudinal ends at x > a/2 (y > b/2) are un- 
der compressive tension in — x(— y)-direction and 
left longitudinal ends at x < —a/2 (y < —b/2) 
are under compressive tension in x (y)-direction 
(see Fig. 1 in Ref. [6]). Note that only in PG 
4 ac'('zz') direction is equivalent to x-direction (y- 
direction). For the lateral ends (for a = they 
are \y\ > b/2 and for a = tt/2 they are \x\ > a/2), 
we used SBC, (only movement of atoms in lateral 
ends in the compression direction is allowed not 
in the ^-direction) and FBC. The FBC (SBC) 
condition is equivalent to suspended graphene at 
two longitudinal ends that is put over a trench 
while it is free at the lateral two ends [10]. Sup- 
ported boundary condition can be created in ex- 
periments by suspending graphene from two lon- 
gitudinal ends and supporting the other two lat- 
eral ends to a substrate which prevents graphene 
to move verticaly [11]. 

We applied a strain rate v = 0.027/ns and 
v = 0.054/ns for a = and a = 7r/2, respec- 
tively. The strain rate is given by v = 2Sx/5000lSt, 
where St = 0.5 fs is our MD-simulation time step, 
Sx = 0.667 pm is used as the compression step after 
each 5000 steps and the factor two is because the 
compression is applied on two longitudinal ends 
(for more details see Ref. [6]). Here I = a or I = b 
if a = or a = 7r/2, respectively. Notice that the 
atoms in the longitudinal ends are fixed during 
each compression step. 

We use the Jarzynski equality [12], i.e. AF = 
— j3~ x ln(exp(— /3W)), which gives a relation between the 
difference of the free energy and the total work done 
on the system (W) during a non-equilibrium evolution 
where f3 = l/fe^T. The averaging is done over several 
realizations of the paths between the initial and the final 
state. We found that averaging over 10 simulations with 
different initial states for each particular case, results in 
a sufficient accurate value for AF (more technical details 
can be found in [7|, [13| ) . 

Elasticity theory predicts that the shape of the lowest 
mode of the buckled state (of a simple bar with length 
/, under axial symmetric load applied at its longitudi- 
nal ends and free from lateral ends), is half a sine wave, 
i.e. Sw = w sin(7rx/l) where Sw is the transverse de- 
flection p, LL4(. For a rect ang ular plate subjected to the 
SBC, elasticity theory [Q, |l4[ predicts the following pos- 
sible deformations 



Sw 



£ 



Wmn sm(n7Tx/a) sin(ra7n//6), 



(i) 



where (m, n) are integers in order to satisfy the SBC and 
Wm n is the amplitude of each mode (m,n). Including 
the appropriate strain energy and using Eq. (pQ) , the mini- 
mum buckling boundary stress for the considered systems 
always occurs for m = 1 and various values of n. It is 
equivalent to a single half wave in the lateral direction 





tt/2 



PG LAGBI LAGBII 
FBC FBC FBC 



0.6% 0.34% 

1.07% 3.0% 



0.01% 
0.2% 



PG LAGBI LAGBII 
SBC SBC SBC 



0.7% 1.45% 0.4% 



TABLE I: Estimated buckling strains for different boundary 
conditions and different a with and without grain boundary. 



and various harmonics n in the direction of compression 
(i.e. perpendicular to the grain boundary). 

From our simulations we found that after many com- 
pression steps GNRs starts to buckle, but the shape of 
the deformed GNRs was found to depend strongly on 
the presence of the grain boundary and on the direc- 
tion of applied compression. Figures [2j a, b,c) (FBC) and 
Figs. [2fd,e,f) (SBC) show snap shots of the deformed 
GNRs without (a,d) and with grain boundary (b,c,e,f), 
beyond the buckling threshold, where e = 2.45% and 
a = 0. The strain is calculated using e = 2Sx/l, 
where I = a (I = b) for a = (a = tt/2). From 
Figs. [2fa,b,c) we see that the deformed shape for PG 
is similar to half a sine wave which is much less the 
case for the LAGBI and LAGBII. The deformations in 
Figs. [2fd,e,f) satisfy the condition m=l in Eq. (pQ), while 
in the direction of the applied compression for LAGBI 
and LAGBII the shape of the deformation is different 
from a sine wave which is most clearly seen around the 
grain boundary. 

The buckling threshold, i.e. e&, is measured by find- 
ing the sudden increase in the average quadratic out of 
plane displacement of the GNR atoms ((h 2 )). The vari- 
ation of (h 2 ), averaged over ten simulations for a = 0, 
versus e are shown in Figs. [3ja,b) for FBC and SBC, 
respectively. The vertical arrows indicate the transition 
points to the buckled state. The buckling strains, e&, are 
listed in Table U for various situations. We found that 
graphene with LAGBII grain boundary subjected 
to FBC vibrates quickly so that the relaxed sys- 
tem (before compression) is buckled and thus the 
buckling strain is zero. This is due to the larger 
angle misorientation (32.2°). However, notice that 
before buckling (FBC) (h 2 ) for PG fluctuates which is 
less for LAGBI. Note that a study of the edge recon- 
struction needs an ab-initio approach which is out of the 
scope of the present study [l5[ . 

Changing a varies e^ significantly. In fact for a = tt/2, 
the buckling strain for LAGBI is three times larger than 
for PG indicating a considerable change in the structural 
deformation of graphene when it is subjected to compres- 
sive tension along the grain boundary. As seen from 
Fig. 02(b), for SBC the largest (smallest) buckling 
strain is for LAGBI (LAGBII) and therefore we 
conclude that graphene with LAGBII is thermo- 
dynamically less stable as compared to LAGBI 
and PG. This is also confirmed by our free en- 
ergy calculations. Note that a larger angle a grain 
boundary results in weaker graphene [16]. 

The obtained buckling strains are comparable to 
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FIG. 2: (Color online) Snap shots of buckled graphene with free lateral boundary condition (a,b,c) and supported 
lateral boundary condition (d,e,f) for perfect graphene (a,d), graphene with LAGBI type grain boundary 
(b,e) and graphene with LAGBII type grain boundary (c,f) where e = 2.45%, a — 0. The red curves give the 
z-deviation averaged over the ^/-direction (to improve visualization the z-components of all atoms were scaled 
by a factor of 3 and the edge atoms were excluded). 
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FIG. 3: (Color online) (a) The out of plane displacement versus applied strain for graphene with free lateral 
boundary condition (a) and supported lateral boundary condition (b) with a = 0. Free energy change (bottom 
panels in (c,d)) and its first derivative (top panels in (c,d)) during compression, for GNRs with (LAGBI, 
LAGBII) and without (PG) grain boundary for FBC (c) and SBC (d). 



the one obtained from recent buckling experiments i.e. 
0.7% [13, LL8|. Our theoretical buckling strains are a lit- 
tle smaller than those found in the experiments, which we 
attribute to the presence of a weak van der Waals inter- 
action between the substrate and graphene in the experi- 
ment. Note that strains are more than an order of magni- 
tude smaller than those where fracture occurs in stretch- 
ing simulations and nanoindentation experiments (they 
were in the range 10-30%) [19j. 

The change in the free energy difference when com- 
pressing the GNRs subjected to FBC (SBC) with a = 
is shown in the bottom panel of Fig. [3jc) (Fig. [31(d)). 
Notice that our non-compressed GNRs (in the beginning 
of the simulations) are flat honeycomb lattice structures 
(and during the first equilibration we did not change its 
size) which are not in a thermomechanical equilibrium 
state at finite temperature. Therefore, the free energy 
of this state should be higher than the one of the equi- 



librium state. The free energy of graphene with 
LAGBII increases with e, i.e. there are no min- 
ima in the free energy curves either for SBC or for 
FBC. This is a confirmation of our previous ar- 
gument about the instability of this system when 
it is suspended. 

Notice that the LAGBI system with FBC (SBC) ex- 
hibits a minimum, i.e. equilibrium state corre- 
sponds to the minimum points in the free energy 
curve, for larger strains, e m =0.73%(1.07%) as com- 
pared to e m =0.52%(0.47%) for PG. The reason is that 
the LAGBI system with a flat surface is much farther 
from thermodynamical stability than the flat PG. The 
top panels in Figs. [3](c) and (d) show the first derivative 
of the free energy for FBC and SBC, respectively. Here 
the transition is continuous because of the finite size of 
the simulated GNR. 

As seen from the bottom panel of Fig. [3fd) at the 



minimum point in the free energy curve the rippled 
state has a lower free energy as compared to the ini- 
tial non-compressed GNRs; i.e. AF(LAGBI) =-10.5 eV, 
AF(PG) =-1.5 eV. Therefore, PG (LAGBI) needs less 
(more) compression steps to reach its equilibrium size. 

In summary, we found that deformations of graphene 
nano-ribbons that are subject to in-plane axial boundary 
stresses is different when the graphene sheet contains a 
grain boundary. In the presence of a grain boundary the 



GNR subjected to compression parallel (perpendicular) 
to the grain boundary has a buckling strain that is 
largest, i.e. 3% (lowest i.e. 0.34%) when the lateral 
edges are free. Large angle grain boundaries result 
in smaller buckling strain and into an instability 
when graphene is suspended. 
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